Study of t] c and rj' c decays into vector meson pairs 



O 

(N 
•*-> 
O 

O 

S 1 

Ql,; 
43 . 



> 

m 

d 
o 



X 



Qian Wang 1 , Xiao-Hai Liu 1 , Qiang Zhao 1,2 
i,) Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, P.R. China 
2) Theoretical Physics Center for Science Facilities, CAS, Beijing 100049, China 

(Dated: October 8, 2010) 

The processes of n c (rj' c ) —¥ VV are supposed to be suppressed by the helicity selection rule (HSR) 
but found to be rather important decay modes in experiment for n c ■ We try to distinguish the 
short-distance transitions via the singly disconnected Okubo-Zweig-Iizuka (SOZI) processes and the 
doubly disconnected (DOZI) processes in n c (n' c ) —> VV. It shows that the SOZI processes can 
be related to the n c (n' c ) wavefunctions at the origin. Therefore, a relation of the decay branching 
ratio fraction between these two decay channels can be established. Such a relation is similar to 
that for J/ip and ip' — > VP, where the so-called "p7r puzzle" has been existed for a while. We 
also show that the intermediate charmed meson loop transitions provide an evading mechanism 
for the DOZI processes. This contribution would turn out to be more important in n' c VV. 
As a consequence, it may produce significant deviations from the SOZI-dominant scenario. Future 
experimental measurement of the n' c — > VV by BESIII should be able to clarify the DOZI-evading 
mechanisms. 

PACS numbers: 13.25.Gv, 11.30.Hv 



I. INTRODUCTION 

Among the charmonium states below DD threshold, there is little information on the rf c decays from experiment. 
Even its mass and width still have large uncertainties as listed in the Particle Data Group p]. From a theoretical 
view point, this state contains rich information on the dynamics in the interplay of perturbative and non-perturbative 
QCD, for which there are still a lot of questions unanswered. 

Our motivation of studying rj c and rj' c simultaneously is driven by two puzzling but interesting questions. Firstly, it 
has been observed that r/ c — > VV, where V stands for light vector meson, is one of the most important decay channels 
for rj c with branching ratios at the order of 10 -3 to 10 -2 . However, these decay channels as well as if c —> VV are 
supposed to be highly suppressed by the so-called helicity selection rule (HSR) 0-11] • Such an observation of the 
HSR violation indicates the importance of QCD higher twist contributions or presence of a non-pQCD mechanism 
that violates the HSR. The contradiction between the data and the HSR expectations based on the perturbative 
QCD (pQCD) has drawn some attention, and various attempts have been made to try to understand the underlying 
dynamics [qHI^. In Refs. [3, HU, we studied the effects of charmed hadron loops as a source of long distance 
contribution which violates the HSR. The results indeed suggest that contributions from intermediate charmed hadron 
loops arc significant and can be taken as an evading mechanism of the HSR. In t] c (r)' c ) — > VV it is natural to expect 
that the charmed hadron loops would contribute and might provide an evading mechanism here. Even for higher 
energy processes such as rjb — > J/ipJ/ip, it was also found that charmed meson loops would enhance the decay rate 
significantly [ll[. In Ref. [lq ], possible contributions from charmed meson loops to e + e~ — > J/ipr] c at W = 10.56 GeV 
were also investigated. 

The second reason that we are interested in rj c (ri' c ) exclusive decays is that they are related with the long-standing 
so-called "pir puzzle" in J/tp(ip') — s- VP. The decays of J/tp and tp' into light hadrons are supposed to be via the 
valence cc annihilations into three gluons in pQCD to the leading order at a typical distance of l/m c . In the heavy 
quark limit, i.e. m c is infinitely large, the decay amplitude will be proportional to the wave function at the origin 
■0(0) (ip'(0)). As a result the following relation is expected for the inclusive final states of light hadrons h, 
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where the phase space factors are taken into account by T to t in the heavy quark limit. For the inclusive decay 
and many exclusive decays, the data exhibit consistencies with this relation rather well, which is called "12% rule" . 
However, some exclusive decay modes, such as p-K and K*K + cc, are found largely violating the "12% rule", which 
is known as the "p7r puzzle" problem. 

Many theoretical efforts have been made in order to understand the origin of such a significant deviation from the 
"12% rule" in J/i(j(ip') — > pn and K*K + cc. In recent works [Tt], [HI: we show that the interferences between the 
strong and EM decay amplitudes in both JjipOtb') — > VP are essential for understanding the "pir puzzle". Similar 
ideas had been proposed in the literature [13, [20(. However, if this is the case, one has to clarify and even quantify 
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the mechanism which suppressed the strong transition amplitude in ip' — > VP and then makes it to be compatible 
with the EM amplitudes in some channels. A numerical study of the overall decay channels for J/ip(ip') —> VP indeed 
suggests such a phenomenon [l8|. In Refs. 0, HJ, it was shown that the intermediate charmed meson loops, which 
serves as a long-distance mechanism for both OZI rule and HSR evasions, can significantly suppress the ip' — > VP 
strong transition amplitudes. Such a mechanism will largely alter the branching ratios for J/ip(ip') — > VP, and makes 
the relation in Eq. ([1} unreliable. 

In fact, the intermediate charmed meson loops could be much more general than we would expect in charmonium 
energy re gion . Further studies of the effects from intermediate meson loops in other processes have been reported in 

Refs. oniin^i!- 

The coherent study of the rj c and r\' c VV will be useful for clarifying the role played by the intermediate meson 
loops and EM transition amplitudes in J/ip(ip') —> VP. Since r\ c and rj c are just the spin partners of J/ip and 
ip' respectively, and they may possess the same spatial wave functions in the heavy quark limit, we would expect a 
similar relation as Eq.([T]) to hold between rj c and rj' cl namely, 



BR( V ' C -» h) = BRjjj -> 77) 
BR(r) c -> h) BR( Vc -> 77) 



Vc(0) 



p»7c 
1 tot 

T V'a 
1 tot 



0.52 ~ 1.56, 



(2) 



where we have taken r]' c (0) / rj c (0) = ip' (0)M(0) = 0.64 (261 ] . and the range of the ratio is displayed considering the 
larger uncertainties of the rj' c total width |l| . Alternatively, we can express this relation as 



R 



Vc(0) 



0.41 



(3) 



which avoids the uncertainties with the total width of rj c . 

In exclusive decays of t] c (t]' c ) — ¥ VV, as discussed earlier the transition amplitudes should be suppressed by the HSR. 
However, since the mass of the charm quark is not heavy enough and due to the non- vanishing light quark masses, 
the HSR violation should occur via both singly disconnected Okubo-Zwcig-Iizuka (SOZI) transitions and doubly 
disconnected (DOZI) processes. The latter can be related to long-distance hadronic transitions via intermediate 
charmed meson loops [22|. For the SOZI transitions, the relation of Eq. §5§ will be respected due to their short- 
distance feature, while the long-distance charmed meson loops will not necessarily respect it. In parallel with the 
observations in J/ip(ip') — > VP, if the charmed meson loops play a more significant role in rj' c — > VV than in rj c — > VV, 
they may violate the relation of Eq. ([3]) and lead to observable phenomena similar to the "p7r puzzle" in these two 
decay channels. 

Our focus in this work is to investigate the role played by the intermediate charmed meson loops in r] c (Vc) VV. 
Also, we mention in advance that due to lack of experimental information to constrain the relative strength between 
SOZI and loop transition amplitudes, some of the results have to be qualitative. But they can be examined by the 
forthcoming BESIII high- statistics experiment. 

As follows, in Sec. II, we discuss the parametrization scheme for r] c (ri' c ) — > VV and constraints from the available 
experimental data. The intermediate charmed meson loops are described by an effective Lagrangian approach in Sec. 
Ill, and numerical results are presented. A brief summary is given in Sec. IV. 



II. PARAMETRIZATION FOR rj (rj' c ) -> VV 

An obvious advantage for the exclusive decays of rj c (r]' c ) — > VV is that these transitions, similar to J/ip(ip') — > VP, 
have only one unique Lorentz structure for the VV P couplings. As stressed a number of times before, this will allow a 
parametrization of the effective coupling constant contributed by different mechanisms. This will help a lot especially 
at this moment when the experimental data for rj' c are not available. 

In Ref. Q, a parametrization scheme is proposed to distinguish the SOZI and DOZI processes by the gluon counting 
rule. Due to the above mentioned property with the VVP coupling, one can factorizc out the DOZI evading amplitude, 
of which the ratio to the SOZI amplitude will give an estimate of the order of magnitude of the contributions from 
the DOZI mechanisms. 

Following Ref. Q , the transition amplitudes for r/ c — > VV can be expressed as 

mV gg \ Vc ) = g 2 R 2 (l + r) 
{uu\V gg \ri c ) = sg(l + 2r) 
(oJ<f>\V gg \ Vc ) - g 2 rRV2 
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(K*+K*-\V gg \ Vc ) = g 2 R 

(p + p-\v gg \vc) = gl ■ (4) 

where V gg is the rj c — > gg — > (qq)(qq) potential, and parameter g$ denotes the coupling strength of the SOZI transitions. 
Parameter r is the ratio of the DOZI transition over the SOZI transition. It should be pointed out that the additional 
gluon exchange in DOZI is not necessarily perturbative. Due to contributions from effective mesonic degrees of 
freedom, namely, intermediate meson exchanges, the DOZI transitions can be evaded by rather soft gluon exchanges 
of which the contributions can be parameterized by r. We also introduce the SU(3) flavor breaking parameter R, of 
which its deviation from unity reflects the change of couplings due to the mass difference between u/d and s. The 
amplitudes for other charge combinations of K*K* and pp are implicated. 

A commonly used form factor is adopted in the calculation of the partial decay widths: 

J- 2 (p)=p 2i cxp(-p 2 /8/? 2 ) , (5) 

where p and / are the three momentum and relative angular momentum of the final-state mesons, respectively, in the 
rj c rest frame. We adopt (3 = 0.5 GeV, which is the same as in Refs. [27l - f30l | . Such a form factor will largely account 
for the size effects from the spatial wavefunctions of the initial and final state mesons. 

The present experimental situation should be clarified. In 2005, BESII Collaboration measured the exclusive decay 
branching ratios of r/ c — > VV (3l|, where they find significant differences from the DM2 results [32| in r\ c — > pp. In 
fact, the PDG averaged value for r\ c — > pp branching ratio is strongly affected by this discrepancy. As a comparison of 
the BESII results and the PDG averaged values, we list in Table [I] the fitting results of BESII data (Fit-II in Ref. 0) 
and PDG2008 Q. The fitted parameters are listed in Table HH 



BR (xl(T 3 ) 


BESII [31J 


Fit-BES 


PDG2008 [1] 


Fit-PDG 


PP 


12.5 ±3.7 ±5.1 


10.7 


20.0 ± 7.0 


11.80 


K*K* 


10.4 ±2.6 ±4.3 


13.6 


9.2 ±3.4 


11.76 


u 


2.5 ±0.5 ±0.9 


2.16 


2.7 ±0.9 


2.23 


OJLJ 


< 6.3 


1.67 


< 3.1 


4.53 


LJ(f> 


< 1.3 


0.33 


< 1.7 


0.02 



TABLE I: The branching ratios for g c -> VV. The data are from BES. Fit-BES is obtained by fitting the BESII data [HI for 
r/c -> 4><j>, K*K* and pp, while Fit-PDG are obtained by fitting the PDG2008 data [J. 



Parameters 


Fit-BES 


Fit-PDG 


r 


-0.16 ±0.15 


0.04 ±0.16 


R 


1.02 ±0.23 


0.91 ±0.16 


9Q 


0.35 ± 0.04 


0.36 ± 0.04 


x 2 


0.5 


2.4 



TABLE II: The parameters determined in Fit-BES and Fit-PDG. 

It is interesting to read that the two fitting results are rather consistent with each other in rj c — > pp, K* K* and 4><j). 
Also, the fitted parameters agree well in these two fits except that the x 2 bas a much larger value in Fit-PDG. The 
fitted branching ratio for pp has significant discrepancies with the PDG averaged value. Differences between these 
two sets of data lead to a relatively large uncertainty with parameter r, while parameters go and R are rather stable. 
It suggests that the DOZI-evading mechanism will only account for an order of r 2 ~ ~ 0.1 in ?/ c — > VV, and the 
SOZI process should be dominant. This observation is within our expectation and consistent with what we find in 
J/il>(tp') — > VP, where the DOZI-evading mechanism in J/ip — > VP is negligibly small fl8j . 

The fitting results also show that the present experimental uncertainties are rather large, and the DOZI-evading 
mechanism cannot be well constrained by the data. It should be pointed out that the relation of Eq. ([3]) may help 
estimate the SOZI contributions in rj' c — > VV given the dominance of the SOZI processes. However, due to lack of data 
for r)' c , it is not possible to estimate the DOZI evading contributions in r]' c — > VV based on the parametrization scheme. 



4 



In order to proceed, we assume that the DOZI-cvading mechanisms are via the intermediate charmed meson loops, 
for which a quantitative study of the effects can be quantified. By taking the upper limit of the DOZI contributions in 
rjc —> VV, we can then estimate the possible impact of the DOZI-evading mechanisms in r)c(r)' c ) — > VV. In particular, 
we would like to examine whether the relation of Eq. §3§ still holds or not with the presence of the DOZI contributions. 



III. DOZI-EVADING MECHANISMS VIA INTERMEDIATE CHARMED MESON LOOPS 



A. Formulation 



We will use an effective Lagrangian approach to estimate the transition amplitudes. In Figs. [TJ [2j [3] and HI the 
Feynman diagrams of rj c decaying to pp, K*K*, wu>, and <p<p via the intermediate charmed meson loops arc presented 
respectively. The relevant Lagrangians based on heavy quark symmetry which describes the coupling between S-wave 
charmonium and the charmed mesons reads [HI, HH , 



C 2 = ig 2 Tr[R c - c H 2ll ^ d ^H u ] + H.c, 
where the S-wave charmonium states are expressed as 



Ran. = 



l-j) 



And the charmed and anti-charmed meson triplet are 



H 2t = -Vm] 



l-j 



(6) 
(7) 

(8) 
(9) 



where V and T>* are pseudoscalar ((D° , D + , D+)) and vector charmed mesons ((D*° , D* + , D* + )), respectively. The 
Lagrangian describing the interactions between light mesons and charmed mesons reads 

C - 



*t/ 



(10) 



with the convention e 0123 = +1, where V and V M denote 3x3 matrices for the pseudoscalar octet and vector nonet 
mesons respectively (35j, 



V 



K+ \ 



V 



7T 

K~ 



V2 + Ve A 



K° 



, V 



V2j V2 

K*° 




(11) 



The following kinematic conventions are adopted, r/ c (p) -> T>(*) (p^V^ (p 3 )[T>(*' ) (p 2 )] -> V(k)V(q), where in 
the square bracket denotes the exchanged charmed meson in the triangle diagrams. There are four types of transition 
amplitudes corresponding to different kinds of charmed meson exchanges in the triangle diagrams. We take the 
amplitudes of r/ c — > (fxp as an example, 



M 



VV'[V] 



d A Pl 
(2tt)4 



(Pi +P2) ■ ekq^e" 



x (pa - P3 ) Q (P1 -P3)x(-g Xp + P M-)\^—F(P$) 



r«3 aia 2 a^ 



(12) 




P 



+ 



P 



iO/ n *0 




+ 




P 



D+(D*+) Vc 



P 




FIG. 1: Feynman diagrams for r\ c — > pp via intermediate charmed meson loops. 
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FIG. 2: Feynman diagrams for r\ c — > K*K* via intermediate charmed meson loops. 
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D + (D*+) 



Vc 



UJ 




FIG. 3: Feynman diagrams for r\ c — ^ uiuj via intermediate charmed meson loops. 



Vc 




FIG. 4: Feynman diagrams for r/ c — > (j><j> via intermediate charmed meson loops. 
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A 2 
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(2tt)* 
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x (-/ A + PlPl 



. P2XP2p-., a 
2 + — 2 ii )(-5« 
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x ("fff + 



/3 



)(-5o 



<A*) 



4^) 



/3 Q 

= l&g^vv fvv've^app Pi{-g x H — )(-.9 P + a - ) 



(15) 



where ai = p\ — m\ , a 2 



Pi 



m%, and 03 



Pi 



The amplitudes of other processes have similar forms as the 



above. We omit them for the sake of brevity. 

Since the couplings in the effective Lagrangians are local ones, ultra-violate divergence in the loop integrals is 
inevitable. We introduce a form factor F(pf) phenomenologically to take into account the non-local effects and cut 
off the divergence in the loop integrals, i.e. 



hpd = n 



A 2 



(16) 



where m, and pi are the mass and four momentum of the corresponding exchanged particle, and the cut-off energy is 
chosen as A; = + ah.QCD with Aqcd = 0.22 GeV [l4|, EH, H^. The value of parameter a is commonly expected 
to be of order of unity. 
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FIG. 5: The partial width fractions T l °°^ vv /r^ 
of a, is also displayed by the dot-dot-dashed line. 



in terms of the parameter a. The ratio of R^^i , which is independent 



B. Numerical results for DOZI-evading contributions 



Before proceeding to the numerical results, we first determine some of the parameters taken in this approach. In 
the chiral and heavy quark limit, the following relations can be obtained [341 l35j. 



9vvv 



gvvv = 



7f ' fvv ' v 



v*v*v 



Affy 
V2 



9v 



my 
m-D* 



ni, lc = 2g 2y / m 11c mv'm'p* 



(17) 
(18) 



where j3 and A are commonly taken as j3 = 0.9, A = 0.56 GeV -1 , while f n is the pion decay constant. 

In principle, the coupling gi should be computed by nonperturbative methods. If we simply estimate the coupling 
t/2 with the vector meson dominance (VMD) argument, it will give gi = y/frhji/ (2mx> /,/,), where and = 405 
MeV being the mass and decay constant of J ftp [33|. This relation gives g Vt ,T>T>* = 7.68, which is a commonly ado pted 
value in the literature. In order to determine g 71 ix>(*>v<-*h we similarly relate it to the coupling of g^'-w = 9.05 [36| . 
which is determined by the experimental data for e + e~ — > DD [37] ]. In Ref. (38j . similar values for S„ e p(»)DW and 
g v "r>(*)T>(*) were adopted. 

The determination of the form factor parameter a will depend on the relative strengths and interferences between 
the SOZI and DOZI-evading amplitudes. As discussed in the previous Section, due to lack of sufficient experimental 
information about the rj' c — > VV, we can only estimate the upper limit of the DOZI-evading contributions based 
on the data for r/ c — > VV. Interestingly, we find that for the same value of the form factor parameter a, the ratio 



r 



loop 

v' c ^vv 



/r" 

/ 7) 



loop 



is rather insensitive to a of a broad range. As mentioned before, the advantage of taking this ratio 



is that the uncertainties of the rf c total width can be avoided. In Fig. [SJ the ratio T °,°^ vv /T °°^ vv in terms of a 
common range of a is displayed. The dot-dot-dashed line (independent of a) denotes the ratio given by Eq. ((3J), which 
turns to be much smaller than the ratio between the charmed meson loops. 
The following points can be learned: 
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r loop {xw- 2 ) 


MeV j) c - 


. vv 


Vc- 


► VV 


U 


0.48 r 


- 0.96 


3.10 r 


- 6.05 


pp 


3.50 r 


- 6.87 


28.04 r 


- 53.70 


K*K* 


3.12 r 


- 6.20 


23.10 r 


- 44.40 




1.15 r 


- 2.28 


9.33 ~ 


- 17.90 



TABLE III: The DOZI-evading contributions to the partial widths of T) c (r)' c ) — > VV estimated by the charmed meson loops. 
The parameter a varies between 0.69 ~ 0.78, which produces about 10% of the experimental data for r\ c — > VV as estimated 
by the parametrization scheme. 



i) In case that the decays of r\ c (r]' c ) — > h to be dominated by the short-distance transitions, i.e. SOZI process, one 
expects that the relation of Eq. would be respected, i.e. 

2 

x T Vc ^ vv ~ 0.41 x T Vc ^ vv , (19) 

where the small DOZI-evading contributions in r\ c — > VV are neglected. 

ii) In case that the DOZI-evading contributions in r/ c — > VV amount to about 10% of the SOZI as shown by the 
fitting results of the previous section, it shows that the intermediate charmed meson loops will be sizeable in rf c — > VV. 
In Table IIII1 we list the DOZI-evading contributions to the partial widths in comparison with the SOZI-dominant 
predictions. These two mechanisms seem to be compatible in r]' c , which may easily violate the relation of Eq. ([3]). Such 
an effect can be directly examined by experimental data for the r)' c VV. Such a possibility that the DOZI-evading 
contributions are sizeable in rf c — > VV exhibits similarities as in tp' — > VP, where the charmed meson loops play an 
important role of causing the deviations from the "12% rule" . 

iii) It should be mentioned again that the decays Tj c (rj' c ) — > VV are strong-interaction-dominant processes, which 
are different from the decays J/ip(ip') — > VP where the EM interaction will also play an important role, especially in 
tp' — > VP [13, E3- We illustrate this point in Fig. [6] By these diagrams, a naive power counting indicates that 

£ for JIM) -> VP, 

~ a e for vM)^W, (20) 

where "T em " and "T str " denote the leading EM and strong transition amplitudes, respectively. It implies that the 
EM contribution in r] c (r]' c ) — > VV will be less important than that in J/i/j(i/j') — > VP. 

An interesting expectation that distinguishes the role played by the EM transitions in r] c (r]' c ) — > VV and J/i/j(ip') — > 
VP is that the long-distance corrections from the charmed meson loops should have the same relative phases in 
different exclusive decay channels. This means that the possible deviations from Eq. (fT9|) caused by the charmed 
meson loops will either enhance or lower the overall rf c — > VV branching ratios. In contrast, in J/i/)(ip') — > VP, 
interferences due to the EM transitions may have different phases in different channels, e.g. the branching ratios of 
ip' — > K*°K° + c.c. is much enhanced in comparison with tp' — > K* + K~ +c.c. [l| . Therefore, observation of systematic 
deviations from Eq. (|19|) would be a strong evidence of contributions from intermediate charmed meson loops. 

IV. SUMMARY 

In this paper, we have investigated the processes of T) c (r)' c ) — > VV. These decay modes are supposed to be highly 
suppressed by the HSR but found to be rather important according to the available data. The intermediate charmed 
meson loop transitions, which are correlated with the DOZI-evading, are introduced to provide a mechanism for the 
HSR violation, and the results indicate that the effect of these loops may be significant, especially in r)' c — > VV. 

Although we are still lack of experimental information to constrain the model parameters in these transitions, we 
show that a parametrization scheme based on the unique VVP coupling structure will allow us to parameterize out 
the SOZI and DOZI-evading transitions in r\ c — > VV. By assuming that the DOZI-evading processes can be recognized 
by the intermediate charmed meson loops as a soft long-distance transition, we find that the SOZI contributions in 
r/c — > VV can be well constrained by the experimental data, and an upper limit for the DOZI-evading contributions 
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Vc(0) 



T 

x em 
Tstr 
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FIG. 6: Schematic diagrams for the EM and strong decays of J/ip(tp') —> VP and rj {ri'^) — > VV. 



can be set. This allows us to proceed and estimate the DOZI-evading contributions in r]' c — > VV. It shows that the 
intermediate charmed meson loops may produce measurable effects in r]' c VV and cause sizeable deviations from 
the relation predicted by the dominance of the SOZI transitions. 

This situation is similar to the exclusive processes of J/ijj(ip') — > VP. There, the intermediate charmed meson 
loops provide a natural explanation for the so-called "pir puzzle". In this sense, the study of rj c (r]' c ) VV would 
provide additional evidence for the intermediate charmed meson loops in charmonium decays, and can be examined 
quantitatively by the future BESIII experiment (39| . 

We should also mention that our conclusions are based on the hypothesis that the flavor components of rj c and r/' c 
are dominated by cc. Thus, the connection of their spatial wavefunctions with those of J/ip and ip' will make sense. If 
r) c or r)' c possesses some other internal structures 0, 0, E3| , the relation between their branching ratios will be affected 
to some extent, which however, is not our focus in this work. 
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